Broadband electromagnetic response and ultrafast dynamics 
of few-layer epitaxial graphene 
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We study the broadband optical conductivity and ultrafast carrier dynamics of epitaxial graphene 
in the few-layer limit. Equilibrium spectra of nominally buffer, monolayer, and multilayer graphene 
exhibit significant terahertz and near-infrared absorption, consistent with a model of intra- and 
interband transitions in a dense Dirac electron plasma. Non-equilibrium terahertz transmission 
changes after photoexcitation are shown to be dominated by excess hole carriers, with a 1.2-ps 
mono-exponential decay that reflects the minority-carrier recombination time. 
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The discovery of graphene - a carbon monolayer and 
building block of graphite, fuUerenes, and nanotubes - 
provides unique opportunities to explore the properties 
of two-dimensional Dirac fermions [TJ. The electromag- 
netic properties and ultrafast carrier dynamics, in par- 
ticular, are important for applications of this new ma- 
terial [21 EJ Hj. In exfoliated graphene, the infrared re- 
sponse is characterized by a universal quantum conduc- 
tivity (jQ = 7re^/2/i, arising from interband transitions 
whose onset energy follows the carrier density [S] [El [TJ [H] . 
A promising route towards large-scale device production 
is epitaxial growth of graphene layers on SiC substrates 
[9]. The optical response of epitaxial graphene is, how- 
ever, much less explored. Recently, first equilibrium 
and time-resolved infrared and terahertz (THz) measure- 
ments were reported on epitaxial graphene with a large 
number (TV = 6-37) of layers [lOl [HI [H] US]. In this 
letter, we present a broadband infrared and ultrafast 
THz study of few-layer epitaxial graphene. Systematic 
thickness variation covers nominally buffer, monolayer, 
and multilayer graphene films. We utilize equilibrium in- 
frared spectroscopy to characterize the broadband con- 
ductivity and transient THz measurements to monitor 
the photoexcited carrier dynamics. This yields momen- 
tum and population relaxation times, and provides in- 
sight into graphene's unusual electrodynamics. 

The samples studied here were grown via thermal de- 
composition on the Si-terminated face of semi-insulating 
6H-SiC(0001) wafers [13]. The film thickness and mor- 
phology was characterized in situ via low-energy electron 
microscopy (LEEM) , showing nanoribbon-like monolayer 
terraces with widths of ~ 60 to 250 nm. Angle-resolved 
photoemission spectroscopy (ARPES) and scanning tun- 
neling microscopy (STM) evidenced a single-crystalline 
character of such domains interspersed with some defects 

[THUS]. 

Measurements of the equilibrium broadband infrared 
response from 20-2500 meV were carried out with 
Fourier- Transform Infrared (FTIR) spectroscopy, using a 
Bruker IFS 66v spectrometer with four different source- 



detector combinations. This was complemented with 
time-domain THz spectroscopy in the low-energy (3-10 
meV) range [17 . The small absorption of the atomically- 
thin layers necessitates optimal suppression of systematic 
errors. We used rectangular SiC wafers (5.5 x 12 mm), 
with graphene growth limited to a central, 4-mm di- 
ameter area via the heater and LEEM cap geometry. 
The sample was then spatially modulated below a 2- 
mm aperture in the spectrometer (inset. Fig. 1), alter- 
nating between graphene (transmission T) and the bare 
substrate (Tq). Systematic errors were thus reduced to 
about ±1% and ±2% for FTIR and THz spectra, respec- 
tively. Figure 1 shows the resulting relative transmission 
spectra T(uj)/Tq{lu) at room temperature, for nominally 
buffer, monolayer, and multilayer graphene. A strong, 
systematic transmission decrease with increasing thick- 
ness is apparent, and all spectra show near-IR absorption 
above ~ 4000 cm^^. Moreover, monolayer and multilayer 
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FIG. 1: (Color online) Relative transmission of nominally 
buffer, monolayer, and multilayer graphene (top to bottom), 
measured by FTIR (solid lines) and time-domain THz spec- 
troscopy (dots). Dotted lines: guide to the eye. Undula- 
tions in the THz data below 100 cm~^ arise from laser drift. 
Hatched: Reststrahlen region of SiC. Inset: sample modula- 
tion scheme. 
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graphene feature far-IR absorption below ~ 500 cm~^, 
with considerable strength given the atomic-scale thick- 
ness. 

More insight is obtained from the optical sheet conduc- 
tivity cr(uj). It relates to the thin- film transmission via 
T{u})/To{uj) + Zoa{u)/{ns + where Zq is the 

vacuum impedance and ng ~ 2.5-3.1 the SiC refractive 
index [TTJ [18] . The influence of the imaginary part of con- 
ductivity is negligible for the given parameters. Thus, we 
can directly obtain the real part ai (w) from the above ex- 
pression, taking into account a frequency-dependent ns 
|19| . Figure 2 shows the resulting sheet conductivity, 
normalized by ctq. The buffer layer response is insu- 
lating (Fig. 2a), as expected from a lack of Fermi- level 
electronic bands |9], with some absorption around 1 eV. 
Figures 2(b) and 2(c) show the mono- and multilayer re- 
sponse after subtracting the buffer layer. It exhibits sig- 
nificant broadband conductivity, with a strong frequency 
dependence in multilayer graphene. 

For further analysis, we calculate ai{u;) for an n-doped 
graphene layer, which at temperature T reads 
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where t is the momentum scattering time and Ep the 
electron Fermi energy [3] [20] ■ The first term is the Drude- 
like intraband conductivity, while the second arises from 
interband transitions above ~ 2iJ|,. For high doping i.e. 
Ep ^ iksT, the conductivity reduces to 
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Here, Ep ~ hv{TrNy/^ , where N is the electron density 
and V ~ 10^ m/s is the Dirac fermion velocity. Thus, 
unhke ordinary conductors, the Drude spectral weight is 
a nonlinear function of N. 

The above model, scaled by the number of layers hl, 
is shown as solid lines in Figs. 2(b) and 2(c). It provides 
for a faithful representation of the measured sheet con- 
ductivity, with riL = 1.5, Ep = 0.45 eV, and r = 2 fs 
for the nominally monolayer sample, and ul — 4.5, 
Ep = 0.22 eV, r = 9.5 fs for multilayer graphene. This 
model comparison shows a consistent scaling of the in- 
traband spectral weight with the high-frequency, inter- 
band response set by n^CQ. The Fermi energy reflects 
the large substrate-induced doping of few-layer graphene, 
corresponding to N — 0.4-1.5 x 10^'^ cm~^, which com- 
pares well to ARPES and STM studies [HIIIl]. We also 
calculated the conductivity with a model that includes 
the gap 2 A seen in ARPES dU [TS]. It is shown by the 
dashed lines in Figs. 2(b) and 2(c), with 2A=250 meV 
and 50 meV respectively. Clearly, at these doping levels 
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FIG. 2: (Color online) Real part of the optical sheet con- 
ductivity (gray), normalized by ctq: (a) buffer layer, (b-c) 
monolayer and multilayer graphene, with the buffer conduc- 
tivity subtracted. Solid and dashed lines: model as explained 
in the text, (d)-(f) 2.5 x 2.5 /xm^ LEEM images of nominally 
buffer, monolayer, and multilayer graphene. 



the effect of the gap on cri(w) is negligible. Figures 2(d)- 
(f) show LEEM images of all three samples, which un- 
derscore their nanoscale inhomogeneity [TS]. Analysis of 
these images yields an average thickness of, respectively, 
0.25, 1.1, and 3.3 ML. This agrees reasonably with the 
fits to the infrared spectra, when considering the lim- 
ited LEEM field of view. The short scattering times 
from the model correspond to a mean-free path VfT = 2- 
10 nm, well below the ~ 100 nm graphene terrace size but 
consistent with scattering from impurities and inherent 
nanoscale ripples in graphene jl6j . 

We now turn to the non-equilibrium THz dynam- 
ics measured via optical-pump THz-probe spectroscopy. 
The graphene layers are excited at room tempera- 
ture with 1.53 eV femtosecond pulses from a 250-kIIz 
Tiisapphire regenerative amplifier and probed with ps 
THz pulses detected via electro-optic sampling [17 . Fig- 
ure 3(a) shows a typical electric field trace Eo{t) trans- 
mitted through the unexcited sample, and the pump- 
induced change AE(t) measured at a fixed pump-probe 
delay At. The sign and negligible phase shift of AE{t) 
indicate a transmission decrease, i.e. added THz con- 
ductivity. Also, AE(t) decays in amplitude with increas- 
ing At but retains its shape (not shown), allowing us 
to determine the overall field change at a fixed time- 
point (arrow. Fig. 3a). The relative transmission change 
is at/To = 2AE/Eo + {AE/E^f, which is plotted in 
Fig. 3(b). The signals peak within the time resolution 
after excitation, followed by a mono-exponential decay 
within the measurable range. 

The maximum incident fiuence of 0.9 ^J/cm^ corre- 
sponds to photoexcited electron and hole densities no = 
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FIG. 3: (Color online) Non-equilibrium THz response, (a) 
Reference field i5o (gray) and pump-induced change AE 
(black) at At = 0.4 ps. (b) Transmission changes (symbols) 
at 0.9 /iJ/cm^ pump fluence. Lines: model of Eq. 3, for elec- 
tron (dashed) and hole contributions (dotted) at Tc — 300 K, 
and the sum (solid lines) with time-dependent Tc. Inset: non- 
equilibrium state. (c),(d) pump fluence dependence of ampli- 
tude AT and recombination time tr for mono- (squares) and 
multilayer (circles). 



Po — 4.6 X 10^°cni~^layer~^, given the graphene inter- 
band absorption 2ZQ(jQ/{ns -I- 1) ~ 1.3%. After exci- 
tation, the excess carriers thermalize with the existing 
plasma on a fs timescale, forming a hot Fermi distribu- 
tion with a temperature increased by AT above the lat- 
tice temperature T^. We can estimate AT from energy 
conservation, Ue{TL, N) + AQ = Ue{TL + AT,N + n) + 
Uh{TL + AT,p), where C/e,/i is the electron or hole gas in- 
ternal energy and AQ is the absorbed pulse energy. For 
Dirac fermions Ue,h = 4,klT^ /{'Kv'^h'^)F2iE'j:i^ /keT), 
where F2 is the second-order Fermi integral. For the high- 
est fluence, this yields AT = 86 K for the monolayer and 
201 K for the multilayer sample. These are upper limits, 
as part of the energy is shed by phonon emission. The 
dense electron gas thus represents a heat bath, and the 
THz response can reflect changes of both majority and 
minority carriers (inset. Fig. 3b). This scenario is differ- 
ent from THz studies of thick multilayer graphene, dom- 
inated by photoexcited carriers in undoped layers [T3j . 

For quantitative comparison, we calculate the non- 
equilibrium conductivity change 
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for weak excitation n ^ N and transient carrier tem- 
perature Tc <C Ep/2kB- This model is shown as lines 
in Fig. 3(b), assuming exponential decay n = p = 



Po exp(— t/r/j) with th =1.3 ps and 1.2 ps for mono- and 
multilayer graphene. Equation 3 includes two parts: (z) 
the conductivity of photoexcited electrons, which follows 
from Eq. 2 with a non-equilibrium electron Fermi energy 
E^p + AE^p. For weak excitation, AEj, ~ {n/2N)E'p 
which renders the electron contribution linear in n. How- 
ever, the electron response - shown as dashed lines in Fig. 
3(b) ~ is about an order of magnitude too small, and thus 
fails to explain the observed signals, {ii) The conductiv- 
ity of photoexcited holes forms the second term in Eq. 3, 
which we derive from the full intraband expression using 
a hole distribution confined to the valence band [TTJ [2T] . 
It exhibits a generally non-linear dependence on the hole 
Fermi energy Ep{p) and, likewise, on the hole density 
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where Li2 is the dilogarithm. The non-equilibrium hole 
response [dotted lines. Fig. 3(b)] by far supersedes the 
electron contribution and yields a close description of the 
THz transmission change. This stark difference between 
electron and hole conductivity, evident in our highly- 
doped layers, reflects graphene 's unusual sensitivity of 
the Drude spectral weight on the carrier distribution |13j . 

Three more aspects should be noted. First, the above 
evaluation assumed Tc = T^. In contrast, the solid 
lines in Fig. 3(b) show the sum of electron and hole 
signals with a time-dependent carrier temperature Tc — 
Tl + ATexp(— i/Tc), with a coohng time Tc = 1.4 ps as 
per Ref. [T^ Clearly, the influence of cooling on the sig- 
nals is minor. Second, stimulated THz emission due to 
interband transitions at the Dirac point has been pre- 
dicted for photoexcited gapless graphene [5]. The lack 
of this effect in our signals reflects the presence of the 
electronic gap. Finally, as shown by the intensity de- 
pendence in Fig. 3(c), the response is nearly linear in p. 
This is conflrmed by the model, see e.g. Fig. 3(b) where 
the hole response follows the electron curve. The linear 
hole response arises at sufficiently small densities where 
^pip) < —kTc, corresponding to p ^ 4 • 10^°cm~^ at 
300 K. In this limit, p ~ 2k%T^ exp{E'l,/kBTc)/iTTV^h'^) 
which renders the hole contribution in Eq. 3 proportional 
to (2v^h/kBTc) X p. 

The transmission changes in Fig. 3(b) thus repre- 
sent the population dynamics of excess holes. Since 
few-layer epitaxial graphene is highly n-doped, electron- 
hole recombination is dominated by the interaction of 
the pump-induced minority carriers (holes) with the 
paramount electron plasma. This explains both the 
mono-exponential kinetics of the transient THz signals 
and the largely excitation-independent effective recombi- 
nation time Tji ~ 1.2 ps in Fig. 3(d). This value of th 
is consistent (within a factor of ~ 3) with calculations of 
Auger and phonon-mediated recombination [J. 

In summary, we studied the broadband equilibrium 
conductivity of few-layer epitaxial graphene - consis- 
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tent with intra- and intcrband transitions of a dense 
Dirac electron plasma - and measured the ultrafast 
minority-carrier recombination time via non-equilibrium 
THz transmission changes. Despite a balance of electron 
and hole excitations, the non-equilibrium THz response 
in these highly n-doped layers is shown to be dominated 
by holes, a confirmation of graphene's unusual electrody- 
namics. 
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acknowledges a scholarship of the Rosztoczy Foundation. 



First, and J. A. Stroscio, Science 317, 219 (2007). 
[17] R. A. Kaindl, D. Hagele, M. A. Carnahan, and D. S. 

Chemla, Phys. Rev. B 79, 045320 (2009). 
[18] M. C. Nuss, K. W. Goossen, J. R Gordon, P. M. 

Mankiewich, M. L. O'Malley, and M. Bhushan, J. Appl. 

Phys. 70, 2238 (1991). 
[19] W. G. Spitzer, D. Kleinman, and D. Walsh, Phys. Rev. 

113, 127 (1959). 
[20] T. Stauber, N. M. R. Peres, and A. K. Geim, Phys. Rev. 

B 78, 085432 (2008). 
[21] V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, Phys. 

Rev. Lett. 96, 256802 (2006). 



[1] A. K. Geim and K. S. Novoselov, Nature Mat. 6, 183 
(2007). 

[2] V. Ryzhii, M. Ryzhii, and T. Otsuji, J. Appl. Phys. 101, 

083114 (2007). 

[3] S. A. Mikhailov and K. Zieglcr, Phys. Rev. Lett. 99, 
016803 (2007); S. A. Mikhailov, Microelectron. J. 40, 
712 (2009). 

[4] F. Rana, Phys. Rev. B. 76, 155431 (2007); F. Rana, P. 
A. George, J. H. Strait, J. Dawlaty, S. Shivaraman, Mvs 
Chandrashekhar, and M. G. Spencer, Phys. Rev. B 79, 
115447 (2009). 

[5] R. R. Nair, P. Blake, A. N. Grigorcuko, K. S. Novoselov, 

T. J. Booth, T. Stauber, N. M. R. Peres, and A. K. Geim, 

Science 320, 1308 (2008). 
[6] Z. Q. Li, E. A. Honrikscn, Z. Jiang, Z. Hao, M. C. Martin, 

P. Kim, H. L. Stormer, and D. N. Basov, Nature Phys. 

4, 532 (2008). 

[7] F. Wang, Y. B. Zhang, C. S. Tian, C. Girit, A. Zettl, M. 

Crommie, and Y. R. Shen, Science 320, 206 (2008). 
[8] K. F. Mak, M. Y. Sfeir, Y. Wu, C. H. Lui, J. A. Misewich, 

and T. F. Heinz, Phys. Rev. Lett. 101, 196405 (2008). 
[9] W. A. do Heer, C. Bcrgor, X. S. Wu, P. N. First, E. H. 

Conrad, X. B. Li, T. B. Li, M. Sprinkle, J. Hass, M. 

L. Sadowski, M. Potemski, and G. Martinez, Solid State 

Commun. 143, 92 (2007). 
[10] J. M. Dawlaty, S. Shivaraman, J. Strait, P. George, M. 

Chandrashekhar, F. Rana, and M. G. Spencer, Appl. 

Phys. Lett. 92, 042116 (2008). 
[11] J. M. Dawlaty, S. Shivaraman. ,J. Strait, P. George, M. 

Chandrashekhar, F. Rana, M. G. Spencer, D. Veksler, 

and Y. Chen, Appl. Phys. Lett. 93, 131905 (2008). 
[12] D. Sun, Z. Wu, C. Divin, X. Li, C. Berger, W. A. de 

Heer, P. First, and T. B. Norris, Phys. Rev. Lett. 101, 

157402 (2008). 

[13] P. A. George, J. Strait, J. Dawlaty, S. Shivaraman, M. 
Chandrashekhar, F. Rana, M. G. Spencer, Nano Lett. 8, 
4248 (2008). 

[14] E. Rollings, G.-H. Gweon, S. Y. Zhou, B. S. Mun, J. L. 

McChcsney, B. S. Hussain, A. V. Fedorov, P. N. First, 

W. A. de Heer, and A. Lanzara, J. Phys. Chem. Solids 

67, 2172 (2006). 
[15] S. Y. Zhou, G. H. Gweon, A. V. Fedorov, P. N. First, W. 

A. de Heer, D. H. Lee, F. Guinea, A. H. C. Neto, and A. 

Lanzara, Nature Mat. 6, 770 (2007); D. A. Siegel, S. Y. 

Zhou, F. El Gabaly, A. V. Fedorov, A. K. Schmid, and 

A. Lanzara, Appl. Phys. Lett. 93, 243119 (2008). 
[16] G. M. Rutter, J. N. Grain, N. P. Guisinger, T. Li, P. N. 



